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OhI We report on a search for the first generation heavy neutrino that is an isosinglet 

(-H I under the standard SU{2)i gauge group. The data collected with the L3 detector at 

center-of-mass energies between 130 GeV and 208 GeV are used. The decay channel 



Ne — >■ eW is investigated and no evidence is found for a heavy neutrino, Ng, in a 
mass range between 80 GeV and 205 GeV. Upper limits on the mixing parameter 



5^ \ between the heavy and light neutrino are derived. 
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Introduction 

In the Standard Model of electroweak interactions pQ, neutrinos are the only fundamental 
fermions which do not have a right-handed component that transforms as an isosinglet under 
the SU{2)l gauge group. However, additional heavy isosinglet neutrinos occur in various 
models that attempt to unify the presently known interactions into a single gauge scheme, such 
as Grand Unified Theories or Superstring inspired models [2j. Several extended electroweak 
models, including left-right symmetric and see-saw models jS] also predict the existence of such 
neutrinos. 

Heavy isosinglet neutrinos can couple to the W and Z bosons through their mixing with 
the light neutrinos. Constraints on isosinglet neutrino mixing were set by several experiments 
IIHO]- Heavy neutrinos were searched for in leptonic decays of mesons and in neutrino beam 
experiments jlj, resulting in stringent upper limits on the square of their mixing amplitude to 
ordinary neutrinos, l^/^P, down to 10~^ in the mass region below 3 GeV. LEP experiments [5^ 
set limits on |^£p of the order of 10^^ to 10~^ for the neutrino mass range from 3 GeV up to 
80 GeV, and the L3 experiment derived the first limits on \Ue\'^ for neutrino masses above the 
W mass p. 

The data used in this analysis were collected with the L3 detector [2j at LEP at center-of- 
mass energies, y/s, between 192 GeV and 208 GeV corresponding to an integrated luminosity 
of 450 pb~^, out of which ~115 pb~^ were collected at ^/s = 206.5 GeV and ~8 pb~^ at 
■\/s = 208 GeV. Final results also include our earlier data recorded at i/i = 133 — 189 GeV [Hj. 

Production and decay 

This search is performed under the assumption that one heavy isosinglet neutrino N^ is associ- 
ated with each generation of light neutrinos with the mixing amplitude Ue. Neither the mixing 
between light neutrinos and higher isodoublet states nor the mixing among light neutrinos are 
considered [HI. 

In e"'"e~ annihilation, single production of heavy neutrinos occurs via the mixing between 
the heavy neutrino and its associated isodoublet neutrino, as presented in Figure 1: 

e"'"e~ — *> N£i>£. 

The corresponding heavy neutrino pair production cross section is suppressed with respect 
to the single production cross section by an additional \Ui\'^ factor, which is expected to be be- 
low 0.1 for a heavy neutrino mass, mN, larger than 80 GeV [H]. The single production proceeds 
through s-channel Z exchange for all generations. In addition, the first generation heavy neutri- 
nos, Ne, which couple to electrons, are also produced through t-channel W exchange. Figure 2 
shows that the t-channel contributions to the total production cross section are dominant and 
the production cross section for Ne can be as high as 0.7 pb. The production cross section 
for N^ and N^ is below the sensitivity of LEP and these heavy neutrinos are not considered in 
the following. 

Heavy isosinglet neutrinos decay via the neutral or charged weak currents: 



N, -^ Zz/, or Np -^ eW. 



The decay into the Z boson is suppressed by the hmited phase space for heavy neutrinos 
with masses close to the the W and Z masses. For masses above 150 GeV the branching ratios 
reach the asymptotic values Br(Ne — > eW) = 2/3 and Br(Ne —>■ Zuc) = 1/3 [THj. 

Event simulation 

Using the full differential cross section jTTj, a dedicated Monte Carlo generator is constructed 
to simulate the production and decay of the heavy isosinglet neutrinos. Subsequent hadronic 
fragmentation and decays are simulated by the JETSET Monte Carlo program J2] ■ The effects 
of the finite width of the produced W and Z bosons as well as initial and final state radiation 
are taken into account. This Monte Carlo program is used to generate several samples of signal 
events with heavy neutrino masses ranging from 80 GeV up to the kinematic limit. For the 
simulation of background from Standard Model processes, the following Monte Carlo programs 
are used: KK2f |I3| (e+e" ^ qq(7)), PYTHIA [12] (e+e" -^ Ze+e", ZZ), KORALZ [HI 
(e+e- -^ r+r-(7)), KORALW |^ (e+e" -^ W+W"), PHOJET [H] (e+e" -^ e+e'qq), 
DIAG36 [HI (e+e- -^ e+e~T+r^), and EXCALIBUR ^ for other four-fermion final states. 

The Monte Carlo events are simulated in the L3 detector using the GEANT [19J and 
GHEISHA [20] programs, which take into account the effects of energy loss, multiple scat- 
tering and showering in the materials. Time dependent detector inefficiencies, as monitored 
during the data taking, are also reproduced. 

Event signatures and selection 

The present analysis concentrates on the decay channel Ng — >■ eW with W — > jets. The signature 
of these events is one isolated electron plus hadronic jets. Since there is only one neutrino in 
the final state, it is possible to reconstruct the invariant mass of the heavy neutrino, that will 
manifest itself as a peak in the invariant mass distribution. Moreover, this decay channel has 
the largest branching ratio varying between 68% and 45% depending on the heavy neutrino 
mass. The dominant backgrounds come from W+W" production with one hadronic and one 
leptonic W decay (92%), qq(7) (5%) and ZZ production (2%). 

The electron identification and jet reconstruction procedures follow the criteria described in 
Reference El The event selection requires at least two hadronic jets plus one isolated electron. 
The visible energy must exceed 70 GeV and the number of reconstructed tracks must be greater 
than 6. The polar angle 6 of the missing momentum has to be in the range 25° < 6 < 155°. The 
visible mass of the event, rriyis, is reconstructed and, to improve the resolution, it is rescaled as: 
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'Pu + E. 



VIS 



where p,y is the missing momentum of the event, and Ems is the visible energy. Figure 3 presents 
the distribution of the rescaled invariant mass, rriresc, after the application of the previous cuts. 
Good agreement is found between data and Monte Carlo expectation. This spectrum is divided 
in two regions of rriresci below and above 100 GeV. In the first, "region 1", the heavy neutrino 
mass is close to the W mass and a significant fraction of W's produced in Ne decays are off- 
shell. For rriresc > 100 GeV, "region 2", the W's are produced mostly on-shell. In this case 
a kinematic fit improves the resolution on the mass measurement, the determination of jet 
energies and angles, and the missing momentum direction for both the signal and the W+W" 



background. Four-momentum conservation and the constraint that the invariant mass of the 
hadronic jets is equal to the W mass, are imposed in the fit. In region 1 we select 27 data 
events with 23.6 ± 0.6 events expected from the Standard Model processes. The corresponding 
numbers for region 2 are 794 and 776.2 ±3.5. Figure 4 displays the distribution of the invariant 
mass of the electron and the missing momentum, rngj/, for events in region 2 after the application 
of the kinematic fit. A clear peak coming from the W"'"W~ background is observed at the W 
mass. 

Finally, the W"'"W~ background is reduced by requiring the invariant mass of the electron 
and missing momentum to be outside the W mass region, mgj, < 70 GeV or mey > 90 GeV, 
which rejects 70% of the background events. Figure 5 shows the invariant mass of the events 
accepted after this cut. We observe a good agreement between the data and expected Standard 
Model background: 233 data events pass the selection with 226.5 ± 1.8 events expected from 
the background, out of which 88% are from W"'"W~ production, 9% from qq(7) production and 
3% from ZZ production. 

Results 

As no signal is observed, the 95% confidence level upper limits on the square of the mixing 
amplitude, |t4p, are calculated from the number of the data and background events ^^. In 
region 1, the total number of selected and expected events is used. In region 2, the number of 
events in data and background for a given heavy neutrino mass rriN is defined as the number 
of events with a reconstructed mass in the interval ?tin ± 2cr. The mass resolution a varies 
from 2 to 2.5 GeV over the investigated mass range. The overall selection efficiency for heavy 
neutrino events varies smoothly from 20% up to 45% depending on the values of mN and y/s. 
The systematic uncertainty on the signal selection efficiency is mainly due to the uncertainty 
in the simulation and reconstruction of the heavy neutrinos (~3%), the signal Monte Carlo 
statistics (~3%), and the energy calibration (~2%). It is estimated to be 5% relative and is 
taken into account in the limit calculation by reducing the selection efficiency by 5%. 

Figure 6 shows the measured upper limits on the mixing amplitude |?7eP as a function of 
the heavy neutrino mass, along with the expected limits as calculated from a large number of 
Monte Carlo experiments. These results are obtained using the whole data sample collected by 
L3 at LEP, and improve upon and supersede our previously published results fHj. 
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a) 




b) 




Figure 1: Feynman diagrams showing the production of isosinglet neutrinos via a) s-channel 
and b) t-channel. Here d, denotes e, /x or r for the s-channel production. 
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Figure 2: Total cross section for single production of heavy isosinglet neutrinos, e"'"e -^ N^z/^, 
as a function of the center-of-mass energy [11]. 
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Figure 3: Distribution of the rescaled invariant mass, 



rrir 



of the event. The points are the 



data, collected at ^/s = 192—208 GeV, and the solid histogram is the background Monte Carlo. 
The shaded histogram is the predicted e^e~ -^ z/Ne signal for the heavy neutrino masses of 
90 GeV and 180 GeV with the mixing amplitude \Ue\ = 0.1. Both histograms are normalised 
to the same luminosity as the data. The split of the spectrum into "region 1" and "region 2" 
is described in the text. 
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Figure 4: The invariant mass, me;/, of the isolated electron and missing momentum. The points 
are the data, collected at y/s = 192—208 GeV, and the solid histogram is the background Monte 
Carlo. The shaded histogram is the predicted e"'"e~ — *> z/Ne signal for a 150 GeV heavy neutrino 
with the mixing amplitude \Ue\ = 0.1. For better visibility, the normalization for the signal is 
scaled by a factor of 2. The arrows indicate the accepted range of rrieu outside the 70—90 GeV 
window. 
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Figure 5: Distribution of the visible invariant mass of the event, m„js, after the kinematic 
fit. The points are the data, collected at a/s = 192—208 GeV, and the solid histogram is 
the background Monte Carlo. The shaded histogram is the predicted e^e~ — >■ uNe signal for 
a 180 GeV heavy neutrino with the mixing amplitude \Ue\ = 0.1. For better visibility, the 
normalization for the signal is scaled by a factor of 2. 
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Figure 6: Observed and expected upper limits at the 95% confidence level on the mixing 
amplitude |f/eP as a function of the heavy isosinglet neutrino mass m-^. 
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